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ABSTRACT 

We present Keck spectroscopic observations of z > 6 Lyman-break galaxy (LBG) candidates in 
the Subaru Deep Field (SDF). The candidates were selected as i'-dropout objects down to z' = 27 
AB magnitudes from an ultra-deep SDF z'-band image. With the Keck spectroscopy we identified 
19 LBGs with prominent Lya emission lines at 6 < z < 6.4. The median value of the Lya rest- 
frame equivalent widths (EWs) is ~ 50 A, with four EWs > 100 A. This well-defined spectroscopic 
sample spans a UV-continuum luminosity range of —21.8 < Myv < —19.5 (0.6 ~ 5 L^y) and a Lya 
luminosity range of (0.3 - 3) x 10 43 erg s" 1 (0.3 - 3 ££ ya ). We derive the UV and Lya luminosity 
functions (LFs) from our sample at (z) ~ 6.2 after we correct for sample incompleteness. We find that 
our measurement of the UV LF is consistent with the results of previous studies based on photometric 
LBG samples at 5 < z < 7. Our Lya LF is also generally in agreement with the results of Lya- 
emitter surveys at z ~ 5.7 and 6.6. This study shows that deep spectroscopic observations of LBGs 
can provide unique constraints on both the UV and Lya LFs at z > 6. 
Subject headings: cosmology: observations — galaxies: high-redshift — galaxies: evolution 



1. INTRODUCTION 

The last decade saw great progress in our understand- 
ing of the distant Universe as a number of objects at 
z > 6 were discovered, i ncluding ga laxies (|Iye et al.l 



20061: iVanzella et all 120111 ), quasars dFan et all 120031: 



Willott et al. 201 



and 
2009). 



7-ray bursts l|Haislip et al.l 
They provide key informa- 



20061: iTanvir et al 
tion to study the formation and evolution of the earli- 
est galaxies, supermassive black holes, and massive stars 
when the Universe was less than one billion years old. 
The observations of these objects such as the detection 
of Gu nn-Peterson troughs in quasar spectra (F an et all 
2006), together with the measurements of the polariza- 
tion anisotropies in t he cosmic microwave background 
(Komatsu ct al. 2009), indicate that at z > 6 we are ap- 
proaching the epoch of cosmic reionization, during which 
the intergalactic space became transparent to Hl-ionizing 
UV photons. 

The first z > 6 galaxies were discovered to be 
Lya emitters (L AEs) at z ~ 6 . 56 using the narrow - 
band technique (|Hu et al.l 120021: IKodaira etall 120031) . 
This technique has been an efficie nt way to fin d high - 
rcdshift galaxies since the work of ICowie fc Hul (fl998h 
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and iRhoads et alj (|2000l ). The number of z > 6 LAEs 
has increased to almost one hundred with a high suc- 
cess rate of spectroscopic confirm ation due to the pres- 
ence of strong Lya emission lines dTaniguchi et al.l | 2005 
Ive et al l I2006t (Kashikawa et al l I2006t IHu et al.l I2OIC ; 



Ota et all 12010: lOuchi et al.l 120101: iCassata et all I20TT 
Kashikawa ct al. 120 111) . Now Lya s urveys are able to 



detect LAE candidates at z = 7.7 (|Hibon et all 120101: 
iTilvi et all 12010) , However, the narrow-band technique 
has its own limitations. Narrow-band filters are built to 
use dark atmospheric windows with little sky OH emis- 
sion. Such windows are rare at the red end of the optical 
range (even rarer in the near-IR), and are also very nar- 
row (~ 100 A), res ulting in small survey volumes. The 
dropo ut technique ( Steidcl fc Hamilton! 1 19931 : IGiavaliscol 
2002) does not have these limitations. It has produced a 
substantial nu mber of Lyman-break galaxy (LBG) candi- 
dates at z > 6 (iBunker et al.ll2004HDickinson et aJJl2004t 
lYan fc WindhorstJ 120041 : iBouwens et all 120081 ), Most re- 
cently with the power of the new HST WFC3/IR cam- 
era, LBG candid ates at z > 7 (up to z ~ 10) are being 
routinely found ([Bunker et al.1 I201O jFinkelstein et al 



2010; Ocs ch et all 12010 Wilkins et al.l I201O lYan et al 
20lOlBouwens et al.l2011l : lLorenzoni et al.ll2011l andref- 



erences therein), although they are mostly too faint to be 
spectroscopically confirmed by current facilities. 

With the large sample of LBG candidates, the galaxy 
UV luminosity function (LF) at z > 6 has been estab- 
lished. The general result is that the faint-end slope 
of the LF at z ~ 6 is very steep, and the characteris- 
tic lum inosity dims significantly from lower reds hifts to 
z ~ 6 (jBouwens et al] 120071 : iMcLure et all 120091 ) , Now 
one of the main concerns is that there is a lack of well- 
defined spectroscopic LBG samples to cross-check the 
LF. As mentioned above, spectroscopic identifications 
of z > 6 LBG candidates are extremely difficult unless 
they are bright and have strong Lya emission lines, in 
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which case spectroscopy could still be costly ([Stark et al.l 
[Ml lOno et al.l [20llt IVanzella et all [20ll . Unfortu- 
nately most of the known LBG candidates at z > 6 are 
in the Hubble Ultra Deep Field (HUDF) due to the abun- 
dance of the high-quality deep data, and thus they are 
very faint. Even the HST W FC3 early release science 
data (|Windhorst et al.l 1201 ID cover only 40-50 arcmin 
and bright candidates are rare. With a great depth over 
an effe ctive area of ~ 876 arcm in 2 , the Subaru Deep Field 
(SDF; Kashikawa et al. 2004) provides a unique field to 
search for relatively bright LBGs. 

The SDF project has been very successful in searching 
for high-redshift galaxies. Taking advantage of an 8-m 
telescope and a prime-focus camera with a large field- 
of-view (FOV, 34' x 27'), SDF has an impressive depth 
(27.5 ~ 28.5 AB mag) in five broad bands BVRi'z' 
over a survey area of one FOV. Especially noteworthy 
is the deep observations with three narrow-band filters, 
NB816, NB921, and NB973, corresponding to the detec- 
tion of LAEs at z ~ 5.7, 6.5, and 7, respectively. So 
far SDF has spe ctroscopically identified ~ 100 LAEs 
at z ~ 5.7 (e.g. iShimasaku et al.l 120061 ) and z ~ 6.5 
(e.g. IKashikaw a et al J 120061 12011 ) , and a few LAEs at 
z ~ 7 (|Ive et al.l 120061 : lOta eHaTl2l)10[ ) . SDF has also 
found a sample of bri ght LBG ca ndidates at z > 6 
(|Shimasaku et al.ll2005ft up to z > 7 (jOuchi et al.ll2009D . 
Five strong LAEs from a list of i'-dropout objects (or 
LBG candidates) have alrea dy been spec trosco pically 
confi rmed at 6 < z < 6.4 (jNagao et all 12004 I2005L 
120071 ). Note that the difference between LAEs and LBGs 
is somewhat arbitrary and there is no clear separation 
line between them, so we simply call the galaxies found 
by the narrow-band technique as LAEs and those found 
by the dropout technique as LBGs. A LBG is also a LAE 
if it is identified to have a strong Lya emission line. 

In this paper we present our deep spectroscopic obser- 
vations of z > 6 LBG candidates in SDF using a new, 
ultra deep z'-band image (29 hour integration) for target 
selection. This image allows us to select candidates down 
to z' = 27 mag, roughl y one mag deeper t han the LBGs 
found by INagao et all (pOOi [20051 l2007h . The struc- 
ture of the paper is as follows. Section 2 briefly describes 
our selection criteria and follow-up observations of galaxy 
candidates. Section 3 presents the results of our spectro- 
scopic observations. We derive the UV-continuum and 
Lya LFs in Section 4, and summarize the paper in Sec- 
tion 5. Throughout the paper we use a A-dominatcd 
flat cosmology with Ho = 70 km s _1 Mpc -1 , fi m = 0.3, 
and 17a = 0.7. Al l magnitudes are on an AB system 
(|Oke fc Gunrifl983l) . 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Selection of Galaxy Candidates at z > 6 

We selected z > 6 galaxy candidates using the SDF 
broad-band images. The SDF public data have a depth 
of B = 28.45, V = 27.74, R = 27.80, i' = 27.43, 
NB921 = 26.54, and z' = 26.62 (3a detection for point 
sources), covering an effective area of ~ 876 arcmin 2 . 
INagao et all (|2004Ll2005Ll2007| ) have used the public data 
to find bright LBGs down to z' ~ 26.1. Recently the 
SDF team has obtained a much deeper z'-band image 
with a total inte gration time of ~ 29 hours and a depth 
of ~ 27.5 mag (|Poznanski et al.l l2007t iRichmond et al.l 



2009). Our candidate selection was based on this deep 
z'-band image together with the public data in the other 
four bluer broad bands. 

We used the traditional dropout technique, i.e., our 
candidates are «'-band dropout objects. The basic crite- 
ria are 

z' < 27 and i' - z' > 1.7. (1) 

Th e color cut is more string e nt tha n i' — z' > 1.5 used 
by INagao et al.l pOOl [20051 [2001 . This is to reduce 
the number of contaminants scattered into the selection 
region due to large i'-band photometric errors. To re- 
move foreground contaminants, we required that the can- 
didates are not detected (< 2a) in three broad bands 
BVR. We also rejected possible z ~ 6.56 LAE candi- 
dates which are relatively bright in the NB921 band with 
respect to their z-band photometry(z' — NB921 < 1), 
since thes e candidates were bei ng targeted in another 
program (jKashikawa et al.1 120111 ). We obtained 499 i'- 
dropouts in the whole SDF field. We then visually in- 
spected each candidate, and removed those with any pos- 
sible detections in any of the BVR bands and those that 
were likely spurious detections (e.g. blended with bright 
stars). The z'-band image is deep enough in most cases 
of our selection. In the extreme case of z' = 27 in which 
the color cut determines i' > 28.7, we visually inspected 
the candidates and simply required that the candidates 
should not be detected in the i'-band (in addition to the 
BVR bands). We generated an i' — 28.7 point source 
(almost all known z > 6 SDF galaxies are point sources 
in SDF images) and put it in a number of blank regions 
of the z'-band image. We cut out these regions and mixed 
them with other blank regions that do not have the sim- 
ulated point source. As a result, more than 95% of the 
regions with the %' = 28.7 source were visually identified, 
so our visual inspection is reliable in this case. Finally, 
196 promising candidates survived for our follow-up spec- 
troscopy. 

2.2. Keck/DEIMOS Spectroscopy 

The follow-up spe ctroscopic observat ions were carried 
out with DEIMOS ([Faber et al.l l200l on the Keck II 
telescope on 25-28 April 2009. The typical seeing was 
1". A total of 79 galaxy candidates from the above were 
covered by six masks, but only 73 of them were observed 
due to slit conflict. There were roughly 100 slitlets per 
mask; most of th e slitle ts were assigned to the targets of 
IKashikawa et al.1 ([20111 ) and various secondary targets. 
We used the 830 lines mm" 1 grating with the order block- 
ing filter OG550. The wavelength coverage is roughly 
from 6000 to 10,000 A. With a 1" slit width, the resolving 
power was ~ 3600. The total integration time per mask 
was ~ 3 hours, broken into individual exposures of 20 or 
30 min. We also observed a spectrophotometric standard 
star BD+28d4211 in long-slit mode with the same grat- 
ing and order blocking filter. The data were reduced with 
the DEEP2 DEIMOS data pipeline based on the spec2d 
IDL packag^B- The DEIMOS flexure compensation sys- 
tem (FCS) failed in the beginning of our observing run, 
which caused problems with the data reduction, as flex- 
ure could happen in both spatial and spectral directions. 

9 The analysis pipeline used to reduce the DEIMOS data was 
developed at UC Berkeley with support from NSF grant AST- 
0071048. 
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We shifted the spectral images along the both directions 
manually based on sky emission lines and flat-field im- 
ages before we fed the data to the pipeline. The pipeline 
corrected for any small residual shifts. 

We extracted the spectra of our targets and flux- 
calibrated the spectra using alignment stars. There were 
typically 4-5 bright (16 < R < 17) alignment stars per 
mask used to align masks. They were put in 4" square 
boxes rather than 1" slits. The advantage of using align- 
ment stars for flux calibration is that alignment stars and 
the targets were observed under exactly the same condi- 
tions such as transparency and airmass. For each mask, 
we first measured the spectral response from the spectra 
of standard star BD+28d4211, and extracted the spectra 
of the alignment stars. We then calculated the count-to- 
flux (erg s _1 cm~ 2 A -1 ) ratios by scaling the spectra 
of the alignment stars to their broad-band photometry 
(i' or z'). The ratios among different alignment stars 
on a mask agree within 0.1 mag. We also incorporated 
slit loss (~ 0.24) into the count-to-flux conversion ratios. 
The slit loss was estimated by assuming a slit width of 
1" and a stable Gaussian PSF of 1" (typical seeing). Fi- 
nally the average of conversion ratios was applied to the 
spectra of other objects in this mask. Because the PSF 
is comparable to the slit width, the slit loss varies with 
varying seeing, possible offsets between targets and slits, 
and even the sizes of targets. We did not correct for these 
minor changes. 

3. RESULTS 

Among six mask images, one (the first one dur- 
ing which FCS failed) had more than twice lower slit 
throughput due to the failure of FCS or other unknown 
reasons. We identified one very bright LAE (out of 14 
candidates) in this mask (No. 19 in Table 1 and Figures 
1 and 2, see the following paragraphs), but we will not 
include this one in the analysis of spatial density and LF 
in the next section. In the other five masks, we identified 
18 galaxies (out of 59 candidates) with prominent Lya 
emission lines. Our 3<r detection limit is about 0.7x 10~ 17 
erg s _1 cm -2 at ~ 8500 A, and is estimated as follows. 
The detection limit depends on the shape of the Lya 
lines, as narrower lines (for a g iven flux) are easier to 
identify. IKashikaw a et al.l (|2011l ) generated two compos- 
ite Lya emission lines for their z ~ 5.7 and 6.5 LAEs and 
found that the two lines were very similar. We create a 
Lya emission model image with the shape of the compos- 
ite z ~ 6.5 LAE profile. We then scale this model image 
and put a number of them onto the spectral images. The 
flux limit is determined by detecting these simulated Lya 
emission from the spectral images. Briefly we found one 
galaxy out of 14 candidates in the mask with very low 
throughput and found 18 galaxies out of 59 candidates 
in the other five masks. The remaining candidates do 
not show any continuum emission nor line emission, so 
we were not able to identify them. 

Figure 1 shows the thumbnail images of these galax- 
ies in two broad bands i' and z' and one narrow band 
NB921. They were barely detected in the i'-band images 
and totally invisible in the BVR bands. Most of them 
were also barely detected in the NB921 band. Figure 2 
shows the DEIMOS spectra and the redshifts of the 19 
galaxies. The spectra have been flux calibrated and are 
placed on an absolute flux scale using alignment stars. 



We can clearly see asymmetry in the emission lines of rel- 
atively bright galaxies. This is the indicator of the Lya 
emission line at high redshift due to strong neutral inter- 
galactic medium (IGM) absorption blueward of the line. 
The non-detection in deep BVR images implies that they 
are not likely low-redshift contaminants. In addition, the 
large wavelength coverage rules out the possibility that 
the detected lines are one of the 11 . [O in] A5007. or 
Ha lines. The high resolution of the spectra also ensures 
that they are not [O Ii] A3727 doublets. In a few cases in 
Figure 2 there are some residual sky lines redward of Lya 
that appear like emission lines. They are not the AGN 
feature Nv A1240, because at z ~ 6 the Nv emission 
line is ~ 170 A away. 

Table 1 lists the galaxy coordinates, redshifts, and z'- 
band magnitudes, as well as other properties that will 
be described below. This galaxy sample spans a redshift 
range 6 < z < 6.4 and a magnitude range of 25.1 < 
z' < 27. Redshift for each galaxy is measured from the 
Lya line center using a Gaussian profile to fit the top 
~ 50% of the line from the peak (the rest-frame Lya line 
center is assumed to be 1216 A). Note that the no. 13 
galaxy in our list is no . 2 in the photometric sample of 
IShimasaku et al.1 (|2005| ). 

We measure observed Lya line fluxes by integrating 
the Lya spectra over rest-frame 1215.2 (=1216-0.8) to 
1217.2 (=1216+1.2) A. We do not correct for IGM ab- 
sorption blueward of the line. Table 1 shows the flux 
measurements and the Lya luminosities derived from the 
observed fluxes. Most of the galaxies in our sample have 
the Lya fluxes in a range of (0.7 ~ 2.4) x 10~ 17 erg 
s _1 cm -2 , comparable to the Lya fluxe s in t he SDF 
z ~ 6.5 LAE sample (jKashikawa et al.1 120061 ). The 
strongest Lya emission line (5.8 x 10 -17 erg s _1 cm~ 2 ) in 
our sample is a s bright as the stro n gest LAE s from large 
LAE surveys oflOuchi et al.l (|2009f ). lHu et al.l (|2010l ). and 
IKashikawa et al.1 (|2011| ). indicating that it represents the 
bright end of LAEs at z > 6. Table 1 also includes the 
star formation rates (SFRs) estimated from the Lya lu- 
minosities by 

SFR(Lya) = 9 x 10~ 43 L(Lya) M yr" 1 , (2) 

which is ba sed on the relatio n between SFR and the Ha 
luminosity (|Kennicutd [T998) and the li ne emission ratio 
of Lya to Ha in Case B recombination (|Osterbock et al.l 
119891 ). The derived SFRs are less than 10 M Q yr" 1 for 
all but one galaxy. 

We estimate the Lya rest-frame equivalent widths 
(EWs) using the observed Lya fluxes and z'-band pho- 
tometry. Although the narrow-band NB921 photometry 
consists of pure continuum flux for objects at z < 6.5, our 
galaxies are mostly very faint in this band, preventing us 
from measuring reliable photometry. The broad z' band 
includes emission from both continuum and Lya. To de- 
compose the z'-band photometry we assu me that the UV 
conti nuum slopes (fx oc A' 3 ) are /3 = — 2 dBouwens et al.l 
120091) . We further assume that the continuum blueward 
of Lya is entirely absorbed by IGM. This is reasonable 
as seen from z ~ 6 quasars (|Fan et al.l 12001 . The only 
free parameter is then the continuum level when an ob- 
served spectrum (continuum+Lya) is scaled to match 
the corresponding z'-band photometry. The results of 
the Lya EWs are listed in Column 8 of Table 1 . Column 
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6 shows -M1300, the absolute AB magnitude of continuum 
at rest-frame 1300 A. The EW measurements are rough 
due to uncertainties from the broad-band photometry, 
Lya fluxes, and UV slopes. We vary the assumed slope 
by 0.5 and the typical change on the measured EWs is 
smaller than 15%. For a few galaxies with strong NB921 
detections, we also derive their EWs based on the con- 
tinua from their NB921 photometry. The results mostly 
agree with the above measurements within 25%. Nev- 
ertheless, there is no doubt that most of the galaxies in 
our sample have large EWs. The median value of the 
EWs is 50 A, with 4 galaxies having EWs higher than 
100 A. They ar e on average smaller than thos e in the 
bright sample of lNagao etall (|2004L 120051 120071 ). whose 
Lya EWs are in the range of 90-240 A. 

4. DISCUSSION 
4.1. Spectral Properties 

We know little about the Lya EW distribution of LBGs 
at z ~ 6. Current facilities can only identify z ~ 6 
LBGs with prominent Lya lines, but it is known that 
most LBGs at low redshift do not have strong Lya emis- 
sion. At z = 2 ~ 3, the Lya EW distribution has been 
well determined ba sed on a spectroscopic sample of more 
than 1000 LBGs (IShaplev et al.ll2lM iReddv et al.ll2008t 
iReddv fc Steidelll2009h . A half of the LBGs in this sam- 
ple have Lya absorption lines instead of emission lines, 
and only 17% have Lya EWs greater than 30 A (with 
- 10 % having EWs > 50 A), while almost all the LBGs 
in our sample have EWs > 30 A. Therefore, LBGs in 
our z > 6 sample have much stronger Lya emission 
on average. Statistics of Lya EWs f or LBGs at z > 5 
has b een tentatively i nvestigated ( e.g. IStark et al.ll2010L 
1201 H ). For example, IStark et al.l (|2011|) found that at 
-21.75 < Muv < -20.25 the fractions of LBGs with 
EWs > 25 and > 55 A are roughly 20% and 8%, respec- 
tively, and at —20.25 < Muv < —18.75 the two fractions 
increase rapidly to ~ 55% and ~ 25%. Our sample spans 
a UV luminosity range of —21.8 < Mi 30 o < —19.5, and 
the two fractions in our sample are 25% (15/59) and 14% 
(8/59), broadly consistent with the trend shown in their 
results. 

The Lya strength may also vary with luminosity. Some 
studies have reported the inverse relationship between 
Lya EW and UV luminosity (e.g. IShaplev et al l 120031: 
lAndo et all 120061 IReddv et alJ I2008t IStark et all 120101 1 . 
while others claimed that t he relation is not ob vious or 
there is no such relation (e.g. lNilsson et aLll2009| ). Figure 
3 shows the Lya EW as a function of the continuum lu- 
minosity Mi 300 for our sample. The dashed line demon- 
strates the detection limit for galaxies at z = 6.2. The 
apparent strong correlation between EW and Mi 300 in 
our sample is likely due to the nature of flux-limited sur- 
veys. Nevertheless, our sample has a similar luminosity 
range as the range of the z ~ 3 LBG sample mentioned 
above. In the ne xt section we will use the EW distri- 
bution at z ~ 3 l|Reddv fc Steidel [2009) to correct for 
sample incompleteness. 

4.2. Sample Completeness 

The galaxies presented in this paper provide a well- 
defined flux-limited galaxy sample down to z' = 27. The 



effective area is ~ 340 arcmin 2 (five DEIMOS mask cov- 
erage) and the redshift range is 6 < z < 6.4. In the 
following subsections we will correct for sample incom- 
pleteness, calculate the spatial density of the galaxies, 
and derive the galaxy LF in this redshift range. 

Due to various selection criteria that we applied, our 
sample is incomplete. The sample completeness is com- 
plicated. Here we correct for incompleteness which 
originates from three major steps, source detection, 
galaxy candidate selection, and spectroscopic identifica- 
tion. The first major incompleteness comes from the 
detection of sources in the SDF image. Because the SDF 
field is crowded, any high-redshift galaxies behind fore- 
ground objects are not detected. It is straightforward 
to calculate this incompleteness by measuring the frac- 
tion of the area occupied by bright objects. There is 
another source of incompleteness due to the fact that 
fainter objects are more difficult to detect. We put a 
large number of simulated galaxies (point sources) in the 
SDF z'-band image and detect them using SExtractor 
(|Bertin fc Arnoutdll~99ll in the way used for our galaxy 
candidate detection. We then measure the completeness 
as a function of magnitude z' by counting the fraction of 
detections. The combined completeness is about 75% at 
z' < 25 and drops to - 65 % at z' = 27. 

The second major incompleteness comes from the mag- 
nitude limit [z' < 27) and color cut (i' — z' > 1.7) that 
we applied to the candidate selection. The last major 
incompleteness is from the fact that our sample is biased 
towards LBGs with strong Lya emission, as discussed 
above. We cannot identify LBGs with Lya fluxes below 
our detection limit in the DEIMOS spectra. We use a 
selection function to correct these incompletenesses. 

The selection function is defined as the probability that 
a galaxy with a given magnitude, redshift, and intrin- 
sic spectral energy distribution (SED) meets the criteria 
of our candidate selection and Lya identification. By 
assuming a distribution for the intrinsic SEDs, we cal- 
culate the average selection probability as a function of 
magnitude and redshift. To do this, we first produce a 
lar ge set of galaxy spe c tra fo r a given (Mi 30 o,z), using 
the lBruzual fc Cha riot (2003) stellar population synthe- 
sis model with the Salpeter IMF. The magnitude Mi 30 o 
is the absolute AB magnitude at rest-frame 1300 A, and 
1300 A is chosen to be close to the z'-band effective wave- 
length for z ~ 6 galaxies. The input model parameters 
will be described in detail in the next paragraph. We 
then apply IGM absorption to the model spectra. The 
neutral IGM fraction increases dramatically from z = 5.5 
to z = 6.5, causing complet e Gunn-Peterson troughs in 
some z > 6 quasar spectra ( Fan et al.ll2006f ). It is thus 
critical to predict i' — z' colors. We calcu late I GM absorp- 
tion in the way used by lFan et al.l (|2001h and I Jiang et alJ 
(2008). Finally we measure the apparent magnitudes 
from the model spectra with the SDF filter transmission 
curves. We also incorporate photometric errors into each 
band. It is particularly important for faint objects, as a 
non-negligible fraction of real galaxies may have been 
scattered out of the selection region due to large photo- 
metric errors. The selection probability for this galaxy 
(M1300, z) is then the fraction of model galaxies that meet 
all our criteria. 

The intrinsic SEDs of galaxies depend on their physical 
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properties such as age, metallicity, and dust extinction. 
However, we know very little about physical properties 
of z > 6 galaxies. We determine input parameters for 
the above synthesis models from our HST and Spitzer 
observations of 20 spectroscopically-confirmed LAEs and 
LBGs at 5.6 < z < 7 in SDF (E. Egami et al., in prepa- 
ration). This sample is bright, and the spectroscopic red- 
shifts remove one critical free parameter for SED mod- 
eling. The HST near-IR data provides rest-frame UV 
photometry to decipher the properties of young stellar 
populations, while the Spitzer IRAC data measure the 
amplitude of the Balmer break and constrain the prop- 
erties of mature populations. We find that there is a 
wide variety of SEDs among these galaxies, such as ma- 
ture galaxies with ages > 100 Myr and young galaxies 
with ages ~ 1 Myr. The dust extinction is low to moder- 
ate, consistent with the general trend that higher-redshift 
galaxies ha ye bluer UV continu um slope and lower dust 
extinction (jBouwens et al. 2009). Based on these results, 
our model parameters are set as follows. At a given red- 
shift (from 5.5 to 7.2), we choose to use a mixed grid of 
six ages [1, 2, 5, 10, 40, 100] Myr and three dust extinc- 
tion values E(B-V) = [0.02, 0.1, 0.3]. The metallicity for 
the models with ages of 1 and 2 Myr is 0.005 Z@, for the 
models of with ages of 5 and 10 Myr is 0.02 Z@, and for 
the models of with ages of 40 and 100 Myr is 0.2 Zq. 
The selection function is not sensitive to these physical 
parameters such as age, metallicity, and dust extinction. 
After we generate the gala xy continuum spectra using the 
iBruzual fc Chariot] (|2003| ) model, we add Lya emission 
lines to the spectra. The Lya emission in our galaxies is 
strong, and has significant contribution (~ 0.4 mag) to 
the z'-band photometry. As we discussed in Section 3, 
we know very little about the statistics of Lya strengths 
at z ~ 6. Therefore we assume that the rest-frame EWs 
of Ly a at z ~ 6 have a dist ribution similar to that of 
z ~ 3 (IReddv fc Steidell[2009l) . 

Figure 4 shows the selection function as a function of 
-M1300 and z. The contours in the figure are selection 
probabilities from 5% to 35% with an interval of 5%. 
The sharp decrease of the probability at z ~ 6 is due to 
the color cut of i' — z' > 1.7. The solid circles are the 
locations of the 19 z > 6 galaxies in our sample. Figure 
4 does not include two small constant incompletenesses. 
One is from slit conflict during DEIMOS slit assignment; 
~ 8% candidates were not allocated slits. The other one 
is due to the existence of strong sky OH lines; emission 
lines that happen to be on these OH lines are much more 
difficult to identify. Although in principle this incom- 
pleteness is also a function of magnitude and redshift, 
we find that 5% is a good approximation for our sample 
in the range of 8500-9000 A fe.g. lStark etalEoWt} . Fig- 
ure 4 does not take into account the rejection of z ~ 6.5 
LAEs either, otherwise there will be a 'dip' at z ~ 6.5 in 
the selection function of Figure 4. 

4.3. UV Luminosity Function 

We derive the volume density of the galaxies at z > 6 
using the traditional 1/V a method. The available volume 
for a galaxy with absolute magnitude -M1300 and redshift 
z in a magnitude bin AM and a redshift bin Az is 

V a = f I p(M 1300 ,z)^dzdM, (3) 
JamJaz az 



where p(Mi3oo, z) is the probability function used to cor- 
rect for all the sample incompletenesses described above. 
Then the spatial density and its statist ical uncertainty 
can be written as ([Page fc CarreralfeOOQf ) 

1/2 

Nqal 1 \ qal / . \ 

P=-?7-> ^{p) ~ ~v7 3 (4) 

where N ga i is the number of galaxies in the bin 
(AM, Az). The magnitude and redshift distributions of 
the galaxies in our sample are shown in Figure 4. We 
measure galaxy densities in one redshift bin 5.8 < z < 6.5 
and three magnitude bins [-19.4, -20.4], [-20.4, -21.4], 
and [-21.4, -23]. The lower limits of redshift and mag- 
nitude are chosen to be our detection limit (probability 
in Figure 4 below 5%). The magnitude upper limit 
23 of the last bin is an assumption, comparable to the 
brightest LBG cand idates in the largest LBG sample of 
IMcLure et all (|2009l ). We choose 6.5 as our upper limit 
of redshift because we were not able to identify LBG can- 
didates at z > 6.6 (see subsection 4.5 for discuss) and the 
LAEs at z ~ 6.56 are being targeting in another program 
(see Section 2). 

Figure 5 shows our measurements of the spatial den- 
sities (the filled circles with error bars) in the three 
magnitude bins at (z) ~ 6.2. It also compares with 
the UV LFs at z = 5 , 5.9, and 6.8 from other studies 
(IBouwens et al.l 120071: IMcLure et all 120091: lOuchi et al.l 
120091: iBouwens et all 120111 ). iBouwens et al.l ([20071) ob- 
tained a large sample of LBG candidates from the HUDF 
and other HST datasets. They derived the UV LFs at 
z = 5 and 5.9 (the blue and green dashed lines in Figure 
5) down to a depth of Mjjv ~ —17 and —18, respectively, 
and found very steep faint-end slopes (a ~ —1.7) at these 
redshifts. IMcLure etaLI (|2009f) extended the HST sam- 
ple by including brighter LBG candidates selected from 
the UKIDSS Ultra Deep Survey (UDS). They found very 
similar LFs at z ~ 5 and 6 (the blue and green dotted 
lines in Figure 4). The z = 6.8 LFs in Figure 5 were 
measure d with the data from the new HST WFC3 IR 
camera (jBouwens et al.ll2011l) and with the deep SDF Y- 
band imaging data (|Ouchi et al.1 120091 ) . It is clear that 
our spectroscopic result is in good agreement with the 
trend of the LFs from z — 5.9 to 6.8 from the previous 
studies. Note that k-correction within a small UV wave- 
length rang e (1300 - 1600 A) is negligible because of the 
blue UV slopes. 

We parametrize the galaxy UV LF at z ~ 6 from our 
spectroscopic sample using a Schechter function, 

cf>(M) = O.41nlO0*lO-°- 4 ( M - M *)( Q+1 )cxp(-lO-°- 4 ( M - M *)), 

where <jf is the normalization, M* is the characteristic 
luminosity, and a is the faint-end slope. Because our 
sample is not deep enough to well determine a, and the 
measurements of a at 5 < z < 6 from previous studies are 
quite rob ust and consistent, we fix a to the value of -1.74 
given by IBouwens et al.l (|2007[ ) . We then fit our data to 
the above function. The best fits are 4>* — (1-7 ± 1.2) x 
10~ 3 Mp c 3 and M* = -19.97 ±0.32, consistent with the 
result of IBouwens et alj (|2007D (<j)* = 1.4 x KT 3 Mpc" 3 
and M* = -20.24) and the result of lMcLure et~aTI (|2009f ) 
{(j)* = 1.8 x 10~ 3 Mpc" 3 and M* = -20.04) at z ~ 6. 
The best-fitting M* = —19.97 indicates that our sample 



6 



spans a large luminosity range 0.6 L* ~ 5 L* across the 
characteristic luminosity L* at z ~ 6.2. The result is 
shown in Figure 5. 

Cosmic variance has a minor effect on our results. The 
six masks covered different parts of SDF, and we did not 
exclusively target denser regions of the LBGs candidates. 
We covered 79 candidates in 340 arcmin 2 , while we se- 
lected 196 promising candidates in the whole SDF field 
(~ 876 arcmin 2 ). The candidate surface density in the 
six masks (0.23 = 79/340) is close to the average density 
of the field (0.22 = 196/876). Therefore our sample is 
a representative sample of SDF. We calculate the uncer- 
tainty originating fr om cosmic variance for a ll the SDF 
candidates using the lTrenti fc Stiavellil ([2008) calculator. 
Within a redshift range of 5.8 < z < 6.5 the uncertainty 
is only ~ 20 %, much smaller than statistical errors in 
cj)* and M*. 

Our galaxy sample is a well-defined spectroscopic sam- 
ple of LBGs at z > 6. The comparison in Figure 5 shows 
that our measurement of the UV LF at z ~ 6 are con- 
sistent with the results from deep HST and UKIDSS 
observations. However, our LF largely depends on the 
selection function described in Figure 4. Our tests show 
that the selection function is not sensitive to the physical 
parameters such as age, metallicity, and dust extinction, 
but is sensitive to the distribution of Lya EWs. Be- 
cause the fraction of strong LAEs among LBGs is small, 
we have applied a large correction in the selection func- 
tion. Figure 5 may indicate that the distribution of Lya 
strength at z ~ 6 is similar to that at z ~ 3. Otherwise 
our LF will be quite diffirent from the previous measur- 
ments of the z ~ 6 LFs. Alternatively, the contamina- 
tion rates in previous LBG samples shoul d be small. For 
examp le, the contamination rate in the iBouwens et al.l 
(2007$ sample was claimed to be as small as ~ 3 %, 
and in the IBouwens et all (|2011l ) sample was found to be 
less than 14 %. If their contamination rates are higher, 
the number densities derived from their samples will be 
lower, in contrary to Figure 5. 

4.4. Lya Luminosity Function 

We also derive the Lya LF using the way very simi- 
lar to what we do for the UV LF. The Lya luminositcs 
L(Lya) are listed in Table 1. We first calculate the selec- 
tion function to correct for incompleteness for our galax- 
ies as a samp le of LAEs. Simula t ed ga laxies are gener- 
ated using the lBruzual fc Chariot! (|2003l ) synthesis model 
as discrbied in Section 4.2. The only difference is the as- 
sumption of the Lya EW distribution. The statistics of 
Lya strengths in LAEs at z ~ 6 has been determined 
in several studies. We use the EW probability distribu- 
tion esti mated from a samp le of photometrically selected 
LAEs of lOuchi et "all (|2008f >. Figure 6 shows the LAE se- 
lection function as a function of L(Lya) and z. The con- 
tours in the figure are selection probabilities from 10% 
to 70% with an interval of 10%. The selection function is 
generally similar to the one for LBGs in Figure 4, except 
that the selection probabilities for LAEs are more than 
three times higher than those for LBGs. This is because 
for LAEs we do not need to apply a large correction due 
to the small fraction LAEs among LBGs. 

The volume densities of the LAEs are computed in 
four luminosity (logL) bins [42.35, 42.60], [42.60, 42.85], 
[42.85, 43.10], and [43.10, 43.50]. The lower limit 42.35 



of the first bin is chosen to be our detection limit (proba- 
bility in Figure 6 below 5%). The brightest bin contains 
only one galaxy and the bin size is arbitrary, so we do not 
include this bin in the following analysis. We estimate 
the Lya LF using a Schechter function, 

(f>{\ogL) = lnWcl)*(L/L*) a+1 cxp{-L/L*). (6) 

Most galaxies in our sample cover the faint end of the 
Lya LF, so we are not able to well determine L* . 
Therefore we fix logL* to be the value of 43.0 (e.g. 



Kashikawa et al.l 2006 


; |Hu 


et al.l 


2011: Kashikawa et al.l 


2011) 


when 



above function. The best fits are </>* = (8.2 ± 3.3) x 
10 -5 Mpc~ 3 and a = -1.67 ± 0.54. Figure 7 shows 
our measurements of the densities (filled circles with er- 
ror bars) and the best model fit (solid curve). As we 
mentioned, the selection function is sensitive to the Lya 
EW distribution, especially in the faint end of the sam- 
ple. The filled triangles in Figure 7 represent the vol- 
ume densities if the assumed EW distribution is increased 
by 30 A (to roughly match th e EW distribution for the 
spectroscopic sample shown in lOuchi et al.l ((2008)). The 
density in the faintest bin is increased by a factor of 
two. In Figure 7 we also show the co mparison with th e 
Lya LF measurements at z ~ 6.5 from IHu et al.l (|2010f) , 
iKashikawa et al.l (120111 ). a ndlOuchi et alj (|2010|T The re- 
sults of iKashikawa et all (|201lD and lOuchi et all (|2010f ) 
agree with each othe r in a wi d e lum inosity range cov- 
ered. The result of IHu et all (|201Q[ ) also agrees with 
them in the bright end of the LF, but show a signifi- 
cant disc r epanc y in the faint end, where the density of 
IHu et all (|2010l ) is about three times lower. The rea- 
son is still unclear. Our LF in the faint end is slightly 
hig her than that of IHu et al.l (l2010fl. but is lower than 
the lKashikawa et al.l (|201lD and lOuchi et al.l (|2010t ) LFs. 
Given the large uncertainties our result is generally in 
agreement with these previous LAE surveys. 

4.5. Did We Miss z > 6.6 LBGs? 

The most distant galaxy in our sample is at z ~ 6.4, 
and we did not find any galaxies at z > 6.6. The reason 
is complicated. The deep z'-band image used in this pa- 
per has securely de tected one of th e most distant LAEs 
known at z — 6.96 (|Ive et al.1 120061 ). We successfully re- 
covered this object during our candidate selection, and 
we also have six good z > 6.6 candidates that have been 
observed. If they were galaxies like the z = 6.96 LAE, 
we should have identified them in the DEIMOS mask im- 
ages. This is illustrated in Figures 4 and 6. The detection 
probability of the z = 6.96 LAE is higher than those for 
a half of the LBGs or LAEs in our sample. However, the 
throughput of DEIMOS drops steadily from 9000 A to- 
wards higher wavelength, and in this range the sky back- 
ground is brighter and the OH lines are denser. These 
issues result in significantly lower signal-to-noise ratios 
in the DEIMOS images, but they were not considered in 
Figures 4 and 6. Therefore it is likely that our spectra at 
> 9000 A are not deep enough to identify Lya emisson 
lines. On the other hand, we cannot rule out the possi- 
bility that there are no z > 6.6 LBGs with strong Lya 
emission in the covered region. Due to the increasing 
neutral fraction of IGM at z > 6, the density of LAEs at 
the bright end may decline rapidly towards higher red- 
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shifts (e.g. Kashikawa et al.ll201lUPentericci et al.ll2011h 
ISchenker et al.H201in . 



5. SUMMARY 

In this paper we have reported the discovery of 19 
LBGs from our deep spectroscopic observations of a sam- 
ple of z > 6 LBG candidates in SDF. The candidates 
were selected using the traditional dropout technique 
from an ultra-deep z'-band image. This image, with a 
total integration time of 29 hours, enables us to select 
galaxies down to z' — 27 mag over a wide field. The 
follow-up spectroscopy was made with Keck DEIMOS. 
The 19 LBGs span a redshift range of 6 < z < 6.4 and a 
magnitude range of 25.1 < z' < 27. They have moderate 
Lya emission line strengths compared to those in LAEs 
at similar redshifts. The median value of rest-frame Lya 
EWs is - 50 A, and four LBGs have EWs > 100 A. 

This well-defined spectroscopic LBG sample spans a 
UV luminosity range of 0.6 ~ 5 i{j V across the UV char- 
acteristic luminosity £{j V ant ^ a Lya luminosity range 
of 0.3 ~ 3 L^ ya across the Lya characteristic luminos- 
ity L^ ya . It thus provides unique constraints on both 
the UV and Lya LFs at z > 6. We correct for sample 
incompleteness from source detection, galaxy candidate 
selection, and spectroscopic identification. In particular, 
for the UV LF we assume that the distribution of Lya 
EWs at z ~ 6 is the same as that at z ~ 3. We then 
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TABLE 1 
Properties of the 19 LBGs 



No. 


RA 
(2000) 


Dec 
(2000) 


Rcdshlft 


z' 


-^1300 


/(Lya) 
(10~ 17 erg s _1 cm~ 2 ) 


EW(Lya) 
(A) 


L(Lya) 
(10 42 erg s- 1 ) 


SFR(Lyo) 
(Mm yr _1 ) 


1 


13:25:18.142 


27:32:32.449 


6.240 


26.99 


-19.80 


0.8 


47 


3.5 


3.2 


2 


13:24:55.589 


27:39:20.772 


6.125 


26.69 


-19.63 


1.4 


90 


5.7 


5.2 


3 


13:24:41.333 


27:43:16.601 


6.394 


26.47 


-20.69 


1.3 


35 


5.9 


5.4 


4 


13:24:36.893 


27:43:41.520 


6.343 


26.38 


-20.74 


1.2 


29 


5.2 


4.7 


5 


13:25:11.086 


27:37:47.957 


6.146 


26.29 


-20.51 


1.0 


23 


3.4 


3.1 


6 


13:24:26.030 


27:16:02.993 


5.992 


26.88 


-19.54 


1.2 


81 


4.7 


4.3 


7 


13:24:20.626 


27:16:40.408 


6.267 


26.91 


-19.90 


1.0 


50 


4.1 


3.7 


8 


13:24:05.894 


27:18:37.703 


6.047 


26.82 


-19.74 


0.7 


41 


2.8 


2.6 


9 


13:24:10.769 


27:19:03.918 


6.038 


26.61 


-19.57 


1.7 


120 


7.2 


6.5 


10 


13:24:34.824 


27:14:18.985 


6.233 


26.98 


-19.67 


1.1 


63 


4.1 


3.8 


11 


13:25:21.053 


27:15:24.023 


6.076 


26.34 


-19.95 


1.9 


92 


7.8 


7.1 


12 


13:25:22.190 


27:21:41.069 


6.214 


26.56 


-19.51 


2.4 


170 


9.6 


8.7 


13 


13:25:27.804 


27:28:58.742 


6.226 


25.31 


-21.77 


1.1 


10 


4.6 


4.2 


14 


13:23:38.616 


27:26:15.623 


6.051 


26.73 


-19.83 


0.7 


41 


3.1 


2.8 


15 


13:23:45.758 


27:32:51.342 


6.313 


25.59 


-21.66 


0.9 


10 


4.2 


3.8 


16 


13:23:42.137 


27:33:33.905 


6.038 


26.05 


-20.58 


1.1 


30 


4.5 


4.1 


17 


13:23:50.484 


27:34:15.078 


6.302 


26.72 


-19.80 


1.6 


101 


7.5 


6.8 


18 


13:24:26.107 


27:18:40.450 


6.131 


26.59 


-19.52 


2.0 


136 


7.7 


7.0 


19 


13:25:21.612 


27:42:28.915 


6.164 


25.10 


-21.25 


5.8 


95 


26.5 


24.1 



Note. — The magnitude M1300 is the absolute AB magnitude of continuum at rest-frame 1300 A. EW(Ly«) is the Lya rest-frame 
equivalent width. The Lya properties are not corrected for IGM absorption. 
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i' z' NB921 i' z' NB921 




Fig. 1. — Thumbnail images of the 19 LBGs in two broad bands i' and z' and one narrow band NB921. The image size is 10" X 10". 
North is up and east to the left. These galaxies were barely detected in the i'-band images and totally disappeared in the BV R bands. 
Most of them were also barely detected in the NB921 band. 
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Fig. 2. — Spectra of the 19 LBGs. The Lya emission line is shown in the center of each panel. The spectra have been flux calibrated and 
are placed on the absolute flux scale. 
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Fig. 2. — Continued. 
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Fig. 3. — Lya EW as a function of the continuum luminosity M1300 for our sample. The filled circles represent the 19 galaxies at 
6 < z < 6.4. The dashed line demonstrates our detection limit for galaxies at z = 6.2. The strong relation between EW and M1300 in our 
sample is very likely due to the nature of flux-limited surveys. 
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Redshift 

Fig. 4. — LBG selection function as a function of Afi3oo and z. The contours are selection probabilities from 5% to 35% with an in terval 
of 5%. The filled circles are the locations of the 19 z > 6 galaxies in our sample. The filled star is the z = 6.96 LAE ((Tyc ct al. 200(|). 
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Fig. 5. — UV LF of LBGs at high redshift. The three filled circles are the spatial densities measured from our sample, and the solid curve 
is the best fit to a Schechter function. As comparison, color-coded dashe d and dotted li nes represent UV LFs from previous studies based 
on ph otometric LBG sample s. The blue and green dashed lines: [Bouwcns ct al. ( 20 0W) . Th e blue and green dotted lines: McL ure et al.l 
(2009). The red dashed line: IBouwens etU] (120111) . The red dotted line: IQuchi et al.l (120091) . 
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Fig. 6. — LAE selection function as a function of L(Lyo) and z. The contours are selection probabilities from 10% to 70% with an inte rval 
of 10%. The filled circles are the locations of the 19 z > 6 galaxies in our sample. The filled star is the z = 6.96 LAE (lye et al. 200BT). 
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Fig. 7. — Lya LF at high redshift. The filled circles and triangles are the spatial densities of our sample at z ~ 6.2 for two different Lyo 
EW distributions (sec Sections 4.3 and 4.4). For the purpose of clarity, the triangles have been shifted 0.02 mag along the horizontal axis. 
The solid cu rves are the bes t fits to a Schechter func tion. The dotted, dashed , and dash-dotted lines represent the Lya LFs of LAEs at 
z ~ 6.5 from lHu etUI 1120101) . iKashikawa et all JMiTT) , and lOuchi et al.1 j20Tot) , respectively. 



